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ABSTRACT 

We investigate the vertical structure and elements distribution of neutrino- 
dominated accretion flows around black holes in spherical coordinates with the 
reasonable nuclear statistical equilibrium. According our calculations, heavy nu- 
clei tend to be produced in a thin region near the disk surface, whose mass frac- 
tions are primarily determined by the accretion rate and the vertical distribution 
of temperature and density. In this thin region, we find that 56 Ni is dominant for 
the flow with low accretion rate (e.g., 0.05M Q s _1 ) but 56 Fe is dominant for the 
high counterpart (e.g., 1M Q s" 1 ). The dominant 56 Ni in the special region may 
provide a clue to understand the bumps in the optical light curve of core-collapse 
supernovae. 

Subject headings: accretion, accretion disks - black hole physics - gamma-ray 
burst: general - nuclear reactions, nucleosynthesis, abundances 



(N 

> 

X 



1. Introduction 

The popular model of the central engine in gamma-ray bursts (GRBs) is named neutrino 



dominated accretion flows (NDAFs, see, e.g..|Popham et al 
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The NDAF involves a hyperaccreting spinning stellar black hole with mass accretion rates 
in the range of 0.01 ~ 10M Q s _1 . In general, low (such as 0.05 M Q s _1 ) and high (1 M® s _1 ) 
mass accretion rate correspond to long (T go > 2s) and short (T 90 < 2s) duration GRBs, which 
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stars or a neutron star and a black hole (jEichler et al 



are originated from collapsar ( IWoosleyill993l; iPaczvnski Il998f) and the merger of two neutron 



198a iPaczvhskil 119911 : iNaravan et al. 



19921 ). respectively. The model can give reasonable explanations to the progenitor and ener- 
getics of fireball in GRBs by the neutrino annihilation or magnetohydrodynamic processes. 
The extreme state is a hotbed to produce heavy elements, and the central region of GRBs 
is an ideal location to supply an extreme ly hot and dense state. Actually, nucleosynthesis 
should also be involved in NDAFs model. iLiu et al.l ( 120071 ) studied the radial structure and 
neutrino annihilation luminosity of NDAFs. We assumed that the heaviest nucleus is 4 He, 
which implies that the numerical value of the electron fraction at the radial outer boundary 
is 0.5. ILiu et al.l (120081) studied the vertical structure and luminosity of NDAFs by the above 
assumption. Kawanaka &: Mineshigd (120071 ) assumed that the inflowing gas is composed pri- 
marily of neutron-rich iron group nuclei, and the electron fraction is 0.42, then studied the 
radial structure and st ability of NDAFs for the different mass accretion rate, using a real- 
istic equation of state (ILattimer fc Swestylll99ll ) in order to properly treat the dissociation 
of nuclei. Unfortunately, they did not provide more information on heavy nuclei. To our 
knowledge, the self-consistent vertical distribution of elements without making a limit for 
the electron fraction has not been previously investigated. 

The purpose of this paper is to investigate the element distribution in the vertical 
direction of NDAFs with detailed neutrino physics and precise nuclear statistical equilibrium 
(NSE). The plan is as follows. In Section 2, with the self-similar assumption in the radial 
direction, we numerically s olve the differential equ ations of NDAFs in the vertical direction 
with the proton-rich NSE (jSeitenzahl et al.l 120081 ) and reasonable boundary condition. In 
Section 3, we present the vertical distribution of physical quantities, such as the density, 
temperature and electron fraction, and show the mass fraction of the main elements at the 
various radii for the different mass accretion rates. Discussion and conclusions are made in 
Section 4. 



2. Equations and boundary condition 
2.1. Equations 



In order to facilitate self-similar process, we adopt the Newtonian potential ij) = —GM/r, 
where M is the mass of the central black hole, and suppose the accretion flow is time- 
independent and axisymmetric in spherical coordinates (r, 9, 0), i.e., d/dt = d/d<p = 0. 



The basic equations are composed of continuity and momentum equations (see e.g., 
Narayan fc Yilll995l ; Kue &: Wangll2005l ). The self-similar assumptions are adopted in the 
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radial direction to simplify the hydrodynamic equations, and the vertical ve locity is ignored 
(vh = 0), then we ob tain the vertical hydrodynamic equations as follows (jGu et al. 
Liu et alJlioiOaiboij l. 



2009 



\v 2 + \c 2 + - r 2 Q K 2 = 0, 



= vacate, (2) 

pau 

3 ac 2 
2 riZ K 

where v r and are the radial and azimuthal components of the velocity, the sound speed 
c s is defined as c 2 = p/p, the Keplerian angular velocity is f2x = (GM/r 3 ) 1 / 2 , and a is a 
constant viscosity parameter. 

Furthermore, according the continuity equation, the mass accretion rate can be written 

as 

/■tt/2 

M = -4vrr 2 / pv r sin 9d9, (4) 

where p is the density, and 9 , it/ 2 are the polar angle of the surface and equatorial plane, 
respectively 

The total pressure is expressed as the sum of four pressures: 

P = Pgas + Prad + Pe + Pu, (5) 

where p gas , p ra d, p e , and p v are the gas pressure from nucleons, the radiation pressure of 
photons, the degeneracy pressure of electrons, and the radiation pressure of neutri n os, re - 



spectively. Detailed expressions of the pressure components were given in iLiu et al.l ( 120071 ) . 
Additionally, we assume the polytropic relation in the vertical direction, p = Kp 4 ^ 3 , where 
K is a constant and 7 = 4/3 is the adiabatic index. 



Considering the energetic balance, the energy equation is written as (ILiu et al.ll2010al . 



2012) 



Qwis Qndv Qui (6) 

where Q v is, Qadv, and Q v are the viscous heating, advective cooling and neutrino cooling 
rates per unit area, respectively. We ignore the cooling of photodisintegration of a-particles 
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and other heavier nuclei (ILiu et al.l l2010aJ . 120121 ). The viscous heating rate per unit volume 
<?vis = v W 2 [d( v <t> I ' r ) I '9r} 2 and the advective cooling rate per unit volume q^v = pv r (de/dr — 
(p/ p 2 )dp/dr) (e is the internal energy per unit volume) are expressed, after self-similar 
simplification, as 



9 apv<p 



(7) 



Qadv 



3 (P ~ Pe)v r 



(8) 



where the entropy of degen erate particles is neglec ted. Thus the vertical integration of 
and Qadv are the following (ILiu et al.ll2010al . 120121 ): 



Q 



V1> - 2 I q vis r sin 9d6 

do 



(9) 



Qadv = 2 q^r sin 9d9. 
'e 



(10) 



The cooling due to the neutrino radiation Q u can be defined as (ILee et al.l 120051 ; ILiu et al. 
20121 ) 



Q V = 2J2 q„ k e~ T »*r sin 6d6 

k J °o 



:n) 



where k represents different types of neutrinos and antineutrinos, and q Vk is the sum of the 
cooling rates per unit volume due to the Urea processes, electron-positron pair annihilation, 
nucleon-nucleon bremsstrahlung, and Plasmon decay, hereafter we represent then with ^ 
2, 3, 4), respectively. The optical depth of neutrino r Uk which includes absorption r a ^ k 
and scattering T Sjl>k , is written as 



V k I 'S,V k 1 



where these two optical depth can be defined as 

Jg2 qi rd9 



rr- 

la,v k ~ 



eo 

7ctT 4 



(12) 



(13) 



T, 



S,Vk 



■ J0o 



(14) 
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where <jj and rij (j = 1, 2, 3, 4) are the cross sections and the number density of protons, neu- 



trons , a-particles, and electrons, respectively (e.g., iKohri et al.l 120051 ; IChen &: Beloborodov 
20071 ). I t is noteworthy that different type of neutrinos involve different processes (for de- 
tails, see lLiu et al.ll2007l ). even so, we still show the sum of all the cooling rates for any type 
of neutrinos or antineutrinos as Equation (13) for simplicity. We ignore the absorption and 
scattering reaction by the heavy nuclei in the model. In order to embody the complicated 
microphysics in NDAFs, we assume that the disk is optically t hin to neut r inos with th e mass 
accretion rate is less than 2M s _1 , which was indicated by iLiu et al.l ( 120071 . 120121 ) . The 
electron fraction can be defined as 



n n + n p 



(15) 



where n p and n n are the total number density of protons and that of neutrons, which has to 
satisfy the the condition of electrical neutrality and the chemical potential equilibrium (e.g., 
Liu et all l2~00~7l ) . 



Liu et al.l ( 120071 ) calculated the electron fraction according to the simple NSE equation, 



the condition of electrical neutrality and a bridging formula of electron fraction that is valid 
in both the optically thin and thick regimes. In this paper, we use the strict NSE equations 
(see Section 2.2) to replace the simple one which assumed that the heaviest nuclei is 4 He. 
Meanwhile, the condition of electrical neutrality still holds, which can be written as 

= n e - - n e +, (16) 

m u 

where m u is the mean mass of nucleus, and n e - and n e + are the numb er density of th e 
electron and positron, which can be given by the Fermi-Dirac distribution (ILiu et al.l 120071 ). 



Furthermore, if Z± and N{ are defined as the number of the protons and neutrons of a 
nucleus X iy and its number density (approximately equals the mass fraction) is n,, then the 
electron fraction can also be written as 



F c = ^n i Z i /^n i (Z i + iV i ). 



(17) 



However, we assume that the disk is optically t hin to neutrinos, the previou s bridging formula 
of electron fraction can be simplified as (e.g., lYuanll2005l ; ILiu et al.l 120071 ) 

2fi e - Q 



If 



n n 



n r 



where fi e and Q are the chemical potential and rest-mass energy difference between neutron 
and proton. n n and n p are the number density of the free neutrons and protons, corresponding 
to nucleus with Z = 0, iV = 1 and Z = 1, iV = in Equation (17). 



-6 - 



2.2. Nucleosynthesis 



In our calculations, the NSE equations are required to be applicable for all the varia- 
tions of the electron fraction. NSE established by all nuclear reactions are in the chemical 
equilibrium. ISeitenzahl et al.l (120081 ) studied proton-rich material in a state of NSE, which 
applies to almost all the range of the electron fraction. The complicated and detailed bal- 
ance has been included under the condition of the equilibrium of chemical potential. They 
showed that 56 Ni is favored in NSE under proton-rich conditions (Y c ~ 0.5) being differ- 
ent from the case of domination by the Fe-peak nuclei with the largest binding energy per 
nucleon that have a proton to nucleon ratio close to the prescribed electron fraction. Par- 
ticularly, t he lower limit of the t emperature in the NSE calculation is identified at about 
2 x 10 9 K ( ISeitenzahl et al.l 120081 ). If the temperature is lower than this limit, the NSE 
solutions will not be reliable. Therefore, in our calculations we assume that all nuclear re- 
actions would cease when the temperature is lower than this limit. The reasons are that 
this limit is near the su rface temperature and the components change little under this limit 
(ISeitenzahl et al.l 120081 ) . The NSE code in proton- rich environments can be downloaded 
from http://cococubed.asu.edu/code_pages/nse.shtml In the NSE, the independent 



variables are the density p, temperature T and electron fraction Y e , which are essential in 
the vertical NDAF model according to the description above. 



2.3. Boundary condition 



Naturally, if the gradient of radiation pressure is larger than the gravity, outflows may 
occur, so we c onsidered that it c ould present a boundary condition for the accretion disk. 
Deviated from iLiu et al.l ( I2010ai ) , here a boundary condition is set to be the mechanical 
equilibrium, like solving the principle of Eddington luminosity. Furthermore, we can ignore 
the radial gradient of radial velocity and pressure at the surface of the disk. There are three 
forces to balance, namely gravity, radiation force and centrifugal force in non-inertial refer- 
ence frame. It should be pointed out that NDAFs are extremely optical thick for photons, 
so the radiation pressure of photons at anywhere can be written as 



Prad 



-aT\ 
3 



(19) 



where a is the radiation constant and T is the temperature of gas. The mechanical equilib- 
rium can be written as ( ILiu et al.l 120121 ) 



Prad 1 0=0 O a T 



2GMm, 



■cot# , 



(20) 
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combined with Equation (19), the surface temperature can be derived as ( ILiu et al.ll2012l ) 

,6GMm 



T 



r^cot^n 



(21) 



where m u is the mean mass of a nucleon, o"t is the Thompson scattering cross. Since the 
NSE and boundary condition have been taken into account, the complete equations of our 
model are established. 



3. Numerical results 



Figure 1 shows the variations of the density p, temperature T and electron fraction 1^ 
with 9. The solid and dashed lines represent the solutions at r = 10r g and 100r 9 , where 
r g = 2GM/ c 2 is the Schwarzschild radius, and the thick and thin lines indicate the solutions 
for M = 0.05M Q s _1 and 1M s _1 , respectively. T he variations of the den sity, temperature 



and electron fraction are similar to the solutions in lLiu et al.l ( 12010ai l2012h . We marked the 



lower limit of the temperature ~ 2 x 10 9 K in NSE eqautions in Figure 1(b). M = 0.05M Q s _1 
and M = 1M & s _1 correspond to Y e around 0.49 and 0.47 near the disk surface, respectively. 
The half-opening angles of the disks are also similar to those in lLiu et al.l (120121 ). which have 
the positive correlation with the accretion rate and radius. 

Figure 2 shows the variations of the mass fraction (also approximately equals the number 
density) of the free neutron and proton, and main elements (include 4 He, 52 Cr, 54 Cr, 54 Fe, 
56 Fe, 56 Ni, and 58 Ni, corresponding to the lines with different colors) with 9 at r = 10r g and 
100r 9 for M = 0.05M Q s _1 and 1M s _1 corresponding to (a)-(d). In addition, according 
to Equation (17), the profiles of Y e can be indicated by Figure 2. Furthermore, the open 
angles in different cases correspond to the lower limit of the temperature in NSE code in 
Figure 1(b), which is the limit of open angle in Figure 2. We consider that the nuclear 
reaction will not occur under this temperature limit. 56 Ni dominates at the disk surface for 
M = 0.05M Q s -1 , and 56 Fe dominates for M = 1M & s _1 , corresponding to Y e around 0.49 
and 0.47, respectively. Other heavy nuclei also appear in these cases. The solutions shows 
the proportion of the nuclear matter increases with radius for the same accretion rate. The 
mass fraction of 56 Ni or 56 Fe near the surface increases with radius. In the middle region, 
4 He is dominant for all the accretion rate. The free neutrons and protons are dominant near 
the equatorial plane of the disk in the hot and dense st ate. Most of th e free protons turn 
into the free neutrons due to the Urea process (see, e.g., ILiu et al.l 120071 ). which causes the 
dominant free neutrons and the decrease of electron fraction. In the model, the accretion rate 
determines the density and temperature. Furthermore, the radial and vertical distribution 
of the density and temperature determines the electron fraction and elements distribution. 
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In simple terms, the change of the electron fraction is inversely associated with the accretion 
rate and radius when the free baryons dominate. 



Discussion and conclusions 



4.1. Discussion 



Supernova light curve bumps have been observed in the optical afterglow of some long 



e.g. : 



Galama et al.lll999l ; IWoosley &: Bloom 



durat ion GRBs, which is driven by the decay of 56 Ni 
20061 ). How to produce massive 56 Ni in supernovae accompanied GRBs is a major problem 
remains unsolved. The possi bilities are that 56 Ni is origina t ed from the central engine trans- 



porte d by the outflow (e.g., iMacFadven fe Wooslev 



1999 



MacFadven 



2003 



Surman et al. 



2011 ) or the explosive burning (e.g., Maeda fe Nomota 2003 ~ lMaeda fe Tominaga 2009 ). 



Moreover, the detecti on of Fe Ka X-ray lines can play an important role in understanding th e 
nature of GRBs (e.g. , lLazzati et al]ll999[ iKallman et aDl2003[ IGou et alJlioOa : lButlerll2007l ) . 
The observations of some X-ray afterglows with Beppo-SA X , ASCA, and Chan dra have 
revealed strong Fe Ka emission lines (e.g.. iPiro et al.ll2000l ). iReeves et al.l (120021 ) reported 
on an XMM-Newton observation of the X-ray afterglow of long-duration GRB 011211. The 
X-ray spectrum reveals evidence for emission lines of Magnesium, Silicon, Sulphur, Argon, 
Calcium, and possibly Nickel, arising in enriched material with an outflow velocity of order 
0.1c. Nevertheless, there are no further observations on the metal lines in GRBs since the 
launch of Swift. There are scarcely any observational evidences on Fe Ka lines appearing in 
the X-ray afterglow of short-duration GRBs. Recently, ma ny works have f ocused on how the 
heavy nuclei are produced by the central engine of GRBs. ISurman et al.l ( 1201 ll ) studied the 
nucleosynt hesis, particula r ly for the progenitor of 56 Ni in the hot outflows from GRB accre- 
tion disks. iMetzger et al.l ( 120111 ) suggested that the composition of ultrahigh energy cosmic 
ray s becomes dominated by heavy nucle i at high energies forming GRBs jet or outflows (also 



see lSigl et al.lll995t iHoriuchi et al.ll2012l ). IMetzger! ( 120121 ) investigated the steady-state mod- 
els of accretion disks produced by the tidal disruption of a white dwarf by a neutron star 
or a stellar black hole, which can produce heavier elements via burning the initial material 
of the white dwarf. They discussed the recently discovered subluminous Type I supernovae 
result from those mergers. 

When a massive star collapses to a black hole, powerful supernova occurs. The newborn 
hyperaccreting black hole may power a GRB. The optical light curve bumps of the supernovae 
accompanied with GRBs is driven by the decay of 56 Ni in the outflows coming from the 
central engine. We have described self-consistently how to produce 56 Ni and other elements 
in the central region of GRBs with the NDAF model. Only in the low accretion rate exactly 
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corresponding to long GRBs, 56 Ni dominates near the disk surface. More daringly, if the 
outflow occurs from the disk surface, which consists of 56 Ni and other heavy nuclei, the 
bumps in supernova light curve can be nat urally generated due to 56 Ni decay in the outflow 
from NDAFs. Actually, iLiu et al.l (120121 ) revisited the vertical structure of NDAFs and 
showed that the possible outflow may appear in the outer region of the disk according to 
the calculations of the vertical distribution of the Bernoulli parameter. Furthermore, We 
noticed that the gradient of radiation pressure may be larger than the gravity in the region 
out of the boundary, which has been represented in Section 2.3. Part of the material near 
the surface would become outflow pushed out by the radiation pressure. We also noticed 
that the description of NDAF is similar to that of slim disk. The high-speed outflow may 
be ge nerated near the di s k surface in the two-dimensional simulat ions of supercritical disk 
(e.g., lOhsuga et al.ll2005l ; lOhsugal 120071 ; lOhsuga k, Mineshigd 120111 ). All these suggest that 



the outflow near the disk surface may appear in NDAF model. 

Nickel existed in the outflows from the surface of NDAF with the low mass accretion rate 
may be the most important source of the long-duration GRB 56 Ni production. Analogously, 
iron and other heavy nuclei may also form the outflow injecting into external environment in 
GRBs. The NDAF model with outflows is necessary to be constructed for further theoretical 
explanation of the bumps in the optical light curve of core-collapse supernovae. 



4.2. Conclusions 

In this paper, we revisit the vertical elements distribution of NDAFs around black 
hole in spherical coordinates with detailed neutrino physics and precise nuclear statistical 
equilibrium and discuss the vertical elements distribution in the solutions. The major points 
are as follows: 

1. The heavy elements distribute near the surface of the disk and the free protons and 
neutrons dominate near the equatorial plane. The distribution of the density and 
temperature determine the nuclear matter distribution. 

2. In the region near the surface of NDAFs, Nickel will be dominant for low mass accretion 
rates, whereas iron will be dominant for high accretion rates. The dominant 56 Ni may 
provide a clue to understand the bumps in the optical light curve of core-collapse 
supernovae. 
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Fig. 1. — Variations of the density p, temperature T and electron fraction Y c with at 
r = 10r 9 (solid lines) and 100r g (dashed lines) for M = mM Q s" 1 , where m = 0.05 (thick 
lines) and 1 (thin lines). 
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Fig. 2. — Variations of the mass fraction of the main elements with 9 at r 
for M = O.O5M s" 1 and 1M s" 1 . 



10r g and 100r 9 



